Chapter 3 
Matrices and Eigenvalue Problems 


re rere tee eee 


3,1 Introduction 


In modern mathiemaiics, matrix theory occupies an important place and has applications in almost all 
branches of engineering and physical sciences, Matrices of order m x n form a vector space and they 


define linear transformations which map vector spaces consisting of vectors in IR” or C” into another 
yector space consisting of vectors in IR” or C” under a given set of rules of vector addition and scalar 
multiplication. A matrix does not denote a number and no value can be assigned to it. The usual rules 
of arithmetic operations do not hold for matrices. The rules defining the operations on matrices are 
usually called its algebra. In this chapter, we shall discuss the matrix algebra and its use in solving 
linear system of algebraic equations Ax = b and solving the eigenvalue problem Ax = Ax. 


3.2 Matrices 


An m Xn matrix is an arrangement of mn objects (not necessarily distinct) in m rows and n columns 
in the form 


Gane Bray ede ein 
An) %m2 +: Ginn 


We say that the matrix is of order m Xx n (mm by n). The objects @)), 42, «+5 ee are hese 
of the matrix.-Each element of the matrix can be a real ora complex num * vu & fanaa a 
Agnes ariailes 9y, apy aes or The, Sea ge mah iss : ie uppercase letters 
Column is called its general element. The matrices ea ae ie ite as EPs baonih 
A,B,C f the matrix is understood, ij 

» B,C, ... etc. When the order 0 ; Dag cca 
. Steripats ofa minut are Teal, © ce sae Bie need ee oh types of 

ix, We define the 
Matrix are com re lled a complex matrix. 
plex it is calle 

Matrices, 
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; is called a row yo 
columns Is ca € 
; ‘ and n co 
Row Vector A matrix of order | x n, that is, it has one TOW 


Ctlor 
Or a row m 


atrix of order n and is written as 


soo ell 
[ay @y> --- Ay], OF [a1 4 n 


in which @); (Or a.) is the jth element, 


n is called a Colum) 
. 3 d one colum 
ah Vector Aaihity opunisen 1, that is, it has m rows an 


vector or a-column matrix of order m and is written as 


by by 

by by 
. b) or 5 

Dan b 


in which bi, (or b;) is the jth element. ; 

The number of elements in a row/column vector is called its order. The vectors are usually denoted 
by boldface lower case letters a, b,c, ... etc. Ifa vector has n elements and all its elements are real 
numbers, then it is called an ordered n-tuple in IR", 
numbers, then it is called an ordered n-tuple in C". 


Rectangular matrix 


whereas if one or more elements are complex 


A matrix A of order m x n,m #n is called a rectangular matrix 
Square matrices A matrix A of order m x n in which m = n, that is number of rows is equal to 
the number of colu i i 


n. The elements a,; 
11> 42> ---, Gy, are called the diagonal elements and the li 


of the matrix. 


Null matrix A matrix A of order m Xn in which al 


; | the elements are Zero is called a null matrix 
or a zero matrix and is denoted by 0, 


Diagonal matrix A square matrix A in which all the off-diagonal elements q., ; ¥% j are zero is 
called a diagonal matrix, For example ij> J ; ‘ 


0 


IS a diagonal Matrix of order n,. 


A diagonal matrix is denoted by D. also Written g 
If all the elements of a diagonal matrj 


NX of order pn a 2 Fy]. 
Sty sare Te equ ; 
is called a scalar matrix of order n. “Qual, that is q 


ii ~ & for all i, then the matrix 
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unit [ 0) I T 
ma 1X i 
Trix T an identity matrix OT ¢ 
b eT 7. 


l 
ne dentity matrix is denoted by I, 


equal matrices Two matrices A = (a.) 4 
: ijJmxn and B = 4 
(i) they are of the same order, that . 3 (Pi)oxg are said to be equal, when 


:; : =p,n= 
(ii) their corresponding elements are equal, that j q and 
5) at 1s ai = bi for all i. j. 


submatrix A matrix obtained b itti 
called a submartix of A. As a Si pees some rows and or columns from a given matrix A is 
ention, the given matrix A is also taken as the slivanactix of A 


3.2.1 Matrix Algebra 
The basic operations allowed on matrices are 


(i) multiplication of a matrix by a scalar, 
(ii) addition/subtraction of two matrices, 


(iii) multiplication of two matrices. 


Note that there is no concept of dividing a matrix by a matrix. Therefore, the operation A/B where 
= 4 Ee —— OOOO EEE 
4 and B are matrices is not defined. 


Multiplication of a matrix by a scalar 


Let @ be a scalar (real or complex) and A = (a;;) be a given matrix of order m x n. Then 


B= aA = a(a;) = (aa;;) for all i and /. (3.2) 


The A det of the new matrix B is same as that of the matrix A. 


Addition/subtraction of two matrices | . 
ces of the same order. Then 


Let A = (a,) and B = (b;) be two matri are oe 
| | = [ j 3.3a 
: i) 
| ee tf ; and j. (3.3b) 
and p = (d,) = B= (aj) — (bij) = (aij b;), tor all i and / 
roca € , i ices of the same 
The order of th trix C or, D is the same aS that of the matrices A and B. Matrices 
er of the new matt! ie oan 
order are said to be conformable for addition/subtrac Lose eae mae . 
ble for 4 j> Ma +++? “P 7 
ITA, A ‘ces which are conforma 
| A), ..., A, are p matrices 
then j ‘ as 
ay, aN + OnA2 fea inte oes: 
Sane: 5 i: A.. The order of the matrix is 
atrices Ais 2°” P 


is called a linear combination of the m 
that BEA ee). 2s = ana 
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Properties of the matrix addition and scalar multiplication 


Let A, B, C be the matrices which are conformable for addition and &, B be scalars, 7 hen 
LA+Be=B+aA. (commutative law) 


21 (A+B +CHA+(B4 C) (associative law), 

3 A+O=A (0 is the null matrix of the same order as A). 
4. A+(~ A) = 0, 5. 7 (A + B)= OA + OB. 
6. (4 BYA = OA + BA, 7. & (BA) = OPA. 


8. 1K A®A and Ox Am), 


Multiplication of two matrices 


he product AB of two matrices A and B is defined only when the number of columns in A js equal 
to the number of rows in B, Such matrices are said to be conformable for multiplication. Let A 


~ (aj) 
H ; pas . i 
be an m x n matrix and B (by) be an n * p matrix. Then the product matrix 
Bi Aa) en iy 
5, Ay tee 4), by, by wre (by) bi, 
ipa lb, | b, 
' 21 22 j~ 2/ 2p 
C™= (c,) = AB = : 
Ay yn ise yy | ; ee 
byt Dy2 Pn lobes: 
| Ay) Any tee Ging 
mA~n n Xp 
is a matrix of order m % p, The general clement of the product matrix C is given by 
n 
ey ™ Aybyy + Apbay + ve + AinPy = De ai Dyp (3.5) 
k=| 


In the product AB, B is said to be pre-multiplied by A or A is said to be post-multiplied by B. 
If A is a row matrix of order | % ” and B is a column matrix of order » x J, then AB is a matrix 
of order | % I, that is a single clement, and BA is a matrix of order n x n, 


Remark 
(a) It is possible that for two given matrices A and B, the product matrix AB is defined but the 
product matrix BA may not be defined, For example, if A is a 2 x 3 matrix and B is 4 
4% 4 matrix, then the product matrix AB is defined and is a matrix of order 2 x 4, whereas 
the product matrix BA is not defined, 


(b) If both the product matrices AB and BA are defined, then both the matrices AB and BA are 
aquare matrices, In general AB % BA, Thus, the matrix product is not commutative. 
If AB = BA, then the matrices A and B are said to commute with each other. 

(c) If AB © 0, then it does not always imply that either A = 0 or B = 0. For example, !et 
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ma a) 
a k 4 and B= } 


ab 
0 0 
then AB = 0 0 
0 o| and BA= 
ax+by 0 A BLA 
(d) If AB = AC, it does not alwa 


ys imply that B = 
(e) Define AS=AXA x A ih saten 


eA (K times). The atri k 
integer k is said to be eater os nen, a matrix A such that A’ 


ent. Th = 0 for some positive 

nt. ers - aaa 

‘ nilpotency of the matrix Ao St Wave OF K for which A® = 0 is called the index. 
Ca Ee Sta SS 

(f) If A? 


= A, then A is called an idempotent matrix. 
Properties of matrix multiplication 


ile If A, B, C are matrices of order m xX n,n X Pp and px q respectively, 
(AB)C = A(BC) 


is a matrix of order m x q. 


then 


(associative law) 


. If A is a matrix of order m x n and B, C€ are matrices of order n X p, then 


A (B + C)= AB+ AC (left distributive law). 


. If A, B are matrices of order m x n and C is a matrix of order n X p, then 
(A + B)C = AC+ BC (right distributive law). 

4. If A is a matrix of order m Xn and B is a matrix of order n x p, then 

a(AB) = A(a@B) = (@A)B 


for any scalar a. 


3.2.2 Some Special Matrices 
We now define some special matrices. 


Transpose of a matrix The matrix obtained by interchanging the pomresnOnG ie rows and columns 
of a given matrix A is called the transpose matrix of A and is denoted by A’ or A’, that is, if 


dy, a 5A0 Gy) 
a M 21 m 
ay, ay) er In 
a> M9 Gin 
wep rth} 1 | Sep 
. a, {mp ge then A = . 
a, ay, Gin 
Qy\ G2 Sr 


; ; T T 
. trices A’A and AA 
If A j “x then A? is ann X m matrix. Also, both the product ma 
1s an m X n matrix, then 


are defined, and 


i n Xn square matrix 
ATA = (n X mlm * p) JS ee 
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AAT = (m X n)(n Xm) is an m X m square matrix. 


and 
bel’: 


A column vector b can also be written as [b, 2 -- 


owing results can be easily verified 


The foll 
a column m 


1. The transpose of 
a TOW matrix. 

2. (A)/=A. 

3. (A+B) =A’ + B’, 

4. (AB)! — B’A!. when the matrices A 

A, is defined, then 


a row matrix is atrix and the transpose of a column Matrix ; 
is 


when the matrices A and B are conformable for addition. 


and B are conformable for multiplication. 


If the product A, Ay ..- 
T 
fAyHouc aul =A; Ace AY. 


~” Remark 2 : 
a,,) of order 1 Xn and a column vector 


The product of a row vector a; = (4j) 4i2 --- 
b; = (0; 3; --- 4 )? of order n x 1 is called the dot product or the inner product of the vectors 


. ny 
a; and b,, that is 


n 


b= YL Gin Pej 


‘ij alka) 
k=l 


which is a scalar. In terms of the inner products, the product matrix C in Eq. (3.5) can be written as 


a,'b, a,b» eo a,°b,, 
a,‘b, a,b, See a,-b 

a = = 3.) 
am b, a,,"b> an b, 


Symmetric and skew-symmetric matrices A real square matrix A = (a;) is said to be 


symmetric, if a,;= a; for all i and j, that is A = A’ 


skew-symmetric, if aj =—a; for all i and /, that is A = — Al. 


~ Remark 3 
(a) In a skew-symmetric matrix A = (a;), all its diagonal elements are zero. 
(b) The matrix which is both symmetric and skew-symmetric must be a null matrix. 
T 
e matrix A Bs 


is always skew- i 
er is symymetic. Therefore, a real square matrix A can be written as t 
ymmetric matrix and a skew-symmetric matrix. That is 


A= maa hang 05, 


2 


twl— 


: 20 
if 4° 
priangl B 
re Zero 


Triangular matric i 
es A square i a /O 
wihenevera <7 Thats q matrix A = (q;) is called a lower triangul fe 
J, that is all elements above t : 


matrix if a; = 0, whenever i > ;, he principal diagonal are zero and an uppe 


that is all the elements below the principal diagonal 4 


Scanned with CamScanner 


Matrices j 
ices and Eigenvalue Problems 3.7 


- gate matrix Let A =(a.) be; 
conjus = (;) be a complex matrix. 


e matrix A = (a;) 


4 Let a. 
qhen, th j) 18 called the conjugate : 


denote the complex conjugate of a,.. 

aa matrix of A, 

ermitian and skew-Hermitian matrice 

; Lge ¢s A complex matrix A js “miti x if 

is oer ters ane x A is called an Hermitian matrix if 
in oe ; an matrix if A = — A? A)’ 

Hermitian matrix is denoted by A” or A* Meare as piles: 


Remark 4 


(a) If A is a real matrix, then an Meices 


skew-Hermitian matrix is same as a skew Fee nage ee ae ia 
: -S 


ymmetric matrix. 


(b) In an Hermitian matrix, all the dia 

: : ; gonal elements are real (let a, =x, + iy,; =@,, gi 
Xj + Wy = Xj — Wy; or y, = 0 for all 7), Cty a ON i 
(c) 


In a skew-Hermiti ix, ¢ 

Kone Tes, ee matrix, all the diagonal elements are either 0 or pure imaginary 
ii Xj ; then a. =-— @.. gives : 
} (endings ay = ~ @ ; gives x, + iy, = ~ (x, — iy) or x, = 0 for all j). 


d) For any complex squ i : ait = a : 
(d) y E quare matrix A, the matrix A +A’ is always an Hermitian matrix and the 


. saul ders ay : 
matrix A~ A” is always a skew-Hermitian matrix. Therefore, a complex square matrix A can 
be written as the sum of an Hermitian matrix and a skew-Hermitian matrix, that is 


A= (A+ AT) + (A~ A). 


Nl 


Example 3.1 Let A and B be two symmetric matrices of the same order. Show that the matrix AB 
is symmetric if and only if AB = BA, that is the matrices A and B commute. 


Solution Since the matrices A and B are symmetric, we have 
AT=A and B’=B. 
Let AB be symmetric. Then 
(AB)’= AB, or B/A’=AB, or BA~ AB. 
Now, let AB = BA. Vaking transpose on both sides, we get 
(AB)! = (BA)! = ATB'= AB. 


Hence, the result. 


3.2.3 Determinants Aon 
te a determinant of order n which is denoted by 


s real if the matrix A is real and may be 
f order n is defined as 


x A of order n, We associa 
has a value and this value 1 


omplex. A determinant 0 


With every square matri 
det (A) or |Aj, The determinant 
teal or complex, if the matrix is ¢ 


ay 412 4 1n 

a>, 42 42n 
det(A) = |Al = ; 

Gy Fn2 ~~ aan 
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n n 


=> CD a,M,= dX 4,4; 


jel j=l 


ne i 
=F G@1)” a, M, = > a, A; 
i=l 


(3.7) 


where M,, and 4;, are the minors and cofactors of a;, respectively. 
We give now some important properties of determinants. 


1. If all the elements of a row (or column) are zero then the value of the determinant ; iS zero, 
2, Al =a 


3. If any two rows (or columns) ate interchanged, then the value of the determinant jg 
multiplied by (—1). 

4. If the corresponding elements of two rows (or columns) are EepRoruenal. to each other 
then the value of the determinant is zero. 

5S. If each element of a row (or column) is multiplied by a scalar @ then the value of the 
determinant is multiplied by the scalar a Therefore, if B is a factor of each element of a row 
(or column), then this factor B can be taken out of the determinant. 

Note that when we multiply a matrix by a scalar a, then every element of the matrix js 
multiplied by a. Therefore, |@A| = a" |A| where A is a matrix of order n. 

6. If a non-zero constant multiple of the elements of some row (or column) is added to the 
corresponding elements of some other row (or column), then the value of the determinant 
remains unchanged. 

7. |A + B] + |A| + |B, in general. 


Remark 5 

When the elements of the jth row are multiplied by a non-zero constant k and added to the 
corresponding elements of the ith row, we denote this operation as R; — R; + KR,, where R; is the 
ith row of |A|. The elements of the /th row remain unchanged See the aericnis of the ith row 
get changed. This operation is called an elementary row operation. Similarly, the operation 
C.-C, + kC;, where C; is the ith column of |Aj, is called the elementary column operation. 
Therefore, Shier sismneuteay row (or column) operations, the value of a determinant is unchanged. 


Product of two determinants 
If A and B are two square matrices of the same order, then 


|AB] = 


: AS p j 
Since |A| = |A‘|, we can multiply two determinants in any one of the following ways 
(i) row by row, 
(iii) row by column, 


(ii) column by column, 
(iv) column by row. 
The value of the determinant is same in each case 


’ 
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gank of a matrix 


if wank of 2 matrix A, denoted by r or ;. 
qherefore. the rank of a matrix is the largest value of r, fe 
martix of A whose determinant is not sexs , for w 


tix A of order n, the rank r = 


(A) is the order of the largest non-zero minor of | Al. 


hich there exists at least one r X r 
7 matrix r < min (m, n). For a square 
The rank of a null matrix is zero and 


gub fA Thus, for an m xj 
n if ; 

mat coal | # 0, otherwise r < 

ie the rank oF matrix 1s 0, then it must be a null matrix : 


gxample 3.2 Find all values of w for which rank of the matrix 


fe SI 
A 0 (ame yO 
a 
—6 If -6 1 


js equal to 3. 


Solution Since the matrix A is of order 4, r(A) < 4. Now, r(A) = 3, if |A| = 0 and there is at least 


one submatrix of order 3 whose determinant is not zero. Expanding the determinant through the 
elements of first row, we get 


jal 0 =i 0 
[Al O Sa iat Pa are =i) = | iu) + 11] -6 
Df 6) wt albert 
= pp = 6p + Mp 6 = (w= 1m = 2)(u — 3). 
Setting |A| = 0, we obtain # = 1, 2, 3. For u = 1, 2, 3, the determinant of the leading third order 


submatrix 
pw -l 0 
|A|=|0  -ll =H #0. 
0 Oe 


Hence, r(A) = 3, when pt = | or > or 3. For other values of p, r(A) = 4. 


3.2.4 Inverse of a Square Matrix 


R is called a 
Let A = (a,,) be a square matrix of order n. Then, A is calle 


(i) singular matrix if |A| = 0, 
(ii) non-singular matrix if |A| # 0. 


) <n and non-singular if its rank 


“<j if i kr(A : 
nis singular if its ran ( pee 


Sanam «ogee ee = see A of order 1 1s said to be invertible, if th 
MA) = n. A square non-singular matrix 


; : uch that 
Non-singular square matrix B of order 7 5 


AB=BA=! 


(3.8) 
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of order n. The matrix B 1s called the inverse matrix o¢ 


hat A’ is the inverse of the matrix A) jf A and ” 


where I is an identity matrix 
write B= A! or A= B |. Hence, we say t 


AtA= AAT =I. 
(3,9) 


The inverse, A” of the matrix A is given by 


I 
aPegteee 3 
A = rs adj(A) Big 


where adj(A) = adjoint matrix of A 
= transpose of the matrix of cofactors of A. . 


Remark 6 
(a) (AB).— BAe: 
We have 

(AB)(AB) | =L 


Pre-multiplying both sides first by A“ and then by B | we obtain 
B'A|(AB)(AB) | = B'(A7!A) B(AB)' = B'A™ or (AB) | = B'A™, 
In general, we have (A,A, ... A)! = A; Any ee Gis 
(b) If A and B are non-singular matrices, then AB is also a non-singular matrix. 
(c) If AB = 0 and A is a non-singular matrix, then B must be null matrix, since AB = 0 can be 
pre-multiplied by A. If B is non-singular matrix, then A must be a null matrix, since AB = 0 can 
_ be post-multiplied by B". 
(d) If AB = AC and A is a non-singular matrix, then B = C (see Remark 1(d)). 
(e) (A+ By zAl+ Be in general. 


Properties of inverse martices 


1. If A exists, then it is unique. 

ae UN oe 

2. Tale -1\T we 

: * : (A’)’. (From (AA™)’ = 1" = 1, we get (A')7A’ = I. Hence, the result). 
- Le = di 

: diag (d,, dy, OOF) Ann)s d,, # 0. Then, Dp = diag (1/d,,, dy, oe Udy): 


. The inverse . 
of a non-singular or 4 
upper or lower tri ix i i an uppe! 
lower triangular matrix. angular matrix is respectively 


6 


. The inverse of ; 
a non-singul babs A : : 
7. (Ay ae = ar Symmetric matrix is a symmetric matrix. 
OF any positive integer n. 
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Zr ee] 
sso 3.3 Show that the matrix A= 15 1 0} satisfies the matrix equation AP 6A 
ORE 3 
ie j= 0 where I is an identity matrix of order 3. Hence, find the matrix (i) A | and (ii) A. 
cancion WP Lave 


OMNES Ss Olas = = 0); 


cubstituting in B = A? — 6A? + 11A— I, we get 


3 6 iO) 2b =o —30 pO =i j-* oF 20 
B=|35 —4 —30|- 90 6 -30|+ |55 1 ())|| = || Cs ae 
BOM asa 20 2AP pot a) Til 38 Ona Ouue! 


0 0 0 
=|0 0 0|= 0.- 
0 0 0 


(i) Premultiplying A’ ; 
atad — 6A TA + 1A7A-A =0 


Or ata A? —6A + TE 


4 -1<5) waste Magee Ne ee 


leg fool PIG NO, tie- <0 lee eee 
Se Aes 0: 6 cele aC ataeaee Sine 
em 3. =e a! O19) AU - 
su (2s 61. = 
ii) A-2 2 S| 527° 5 160 
(ii) A? = (A eel 6 os re 55 21 19 


5 asc xls greene 


We can also write 
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3.2.5 Solution of 2 x n Linear System of Equations 
Consider the system of 7 equations in 7” unknowns 
ay\\*| + a\2*2 =n ore hk Ayn Xn = b, 


aX ef 79% an ate AgnXp =F by 


Ay) Xy + a,2X + on. F Ayn Xn ~ by. B11 
In matrix form, we can write the system of equations (3.11) as 
go (3.12) 
a 2 a, b %; 
a, a Ay by Xp 
where A= Repeal =| 
ai G2 athe Qin b,, Xe 


and A, b, x are respectively called the coefficient matrix, the right hand side column vector and the 
solution vector. If b #0, that is, at least one of the elements b,, by, ..., is not zero, then the system 


of equations is called non-homogeneous. If b = 0, then the system of equations is called 
homogeneous. The system of equations is called consistent if it has at least one solution and 


inconsistent if it has no solution. 
Non-homogeneous system of equations. 
The non-homogeneous system of equations Ax = b can be solved by the following methods. 
Matrix method 
Let A be non-singular. Pre-multiplying Ax = b by A7|, we obtain 
x= A7'b. (3.13) 


The system of equations is consistent and has a unique solution. If b = 0, then x = 0 (trivial solution) 
is the only solution. . 


Cramer’s rule 


Let A be non-singular. The Cramer’s rule for the solution of Ax = b is given by 


= |A,| od (3.14) 


where |A.| is the determinant of th i ; righ 
hand side column vector b. e matrix A; obtained by replacing the ith column of A bY the 


We discuss the following cases, 
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stem of equati : 
10NS is ‘ 
Gg Cons F 
4 Eq. (3.14). Stent and the unique solution is obtained by 


sg When |A| = 0 and one or more of jay p_ 4 

? i at j ib 6 = ; 

tions has no solution, that is the system js inconsictony 7 Ste HOt zero, then the system of 
¢ : nt. 


se 3 When |A| = 0 and all lAH=09=4 ; 


: ii Ss sby | 
“4s infinite number of solutions. The System of »N, then the system of equations is consistent 
an 


“quations has at least a one-parameter family of 
pomogencous system of equations 
consider the homogeneous system of equations 


noe (3.15) 
Trivial solution x = 0 is always a solution of this system. 


fA is non-singular, then again x = A™ 0 = 0 is the solution, 


Therefore, a homogeneous system of equations is always consistent. We conclude that non-trivial 


lutions for Ax = 0 exist if and only if A is singular. In this case, the homogeneous system of 
equations has infinite number of solutions. 


Example 3.4 Show that the system of equations 


1 -l 1 || x 4 
2 1 -3ily|=1] 0). 
] 1 || z 2 


has a unique solution. Solve this system using (i) matrix method, (ii) Cramer’s rule. 


Solution We find that 


a Sed 
-l)= 0. 
|Al = |2 | 3) = 1(1 + 3)-2C1-1)+1G 1)=10 # 


(ji4e uw 


T 
solution. Let x = [x, y, 2]. 


Wi Diged ie 
= £3 ty See and b = | 0}. 
(i) We obtain A a A 

1 2 3 


BS es 
Therefore, > 10 3\l2 


Hence, a 2. y = | and z= 1. 
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Tce Sa 
(ii) We have l= 110 ns =P S) 0 + 2(3 — 1) = 20. 
2 I ] 
byt! l 
|A}=|2 0 -3/= 1(0 + 6) — 244 — 2) + 112-0) =_ 49 
il 2 l 
1 -l 4 
|Aj)=|2 1 0) = 1(22—-0)-2- 2-4) + 10-4) = 10, 
Ae bla bes 7A 
|A\| |A.| |A3| 
Therefore, C= — =2, p= ~~ =-1,z= =], 
Jalil “aticanee fay \Al 
Example 3.5 Show that the system of equations 
1 4 3i)x/ [3 
23. Aisa =| 2 
3 2 4]) 2x 5 
has infinite number of solutions. Hence, find the solutions. 
Solutions We find that 
i all 3 B= 3 
Al=Ml2/ese (0, JA\f="2in3 I= 0. 
24 Sane 
i Bo lentes 
Bole, 2a tlis=n0, 0 /An al) 3. o.e 9 
5 4 B= eS 


Therefore, the system of equations has infinite number of solutions. Using the first two equation 


IG Sd) ies 
and solving, we obtain x, = (11 — tion 


10 x,)/5 x Slee oa SISO) 
satisfies the third equation. pass iRosies DO Where e aay: 
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gxample 3.6 
4-9 3x 
2 oral alley 
2 A6t Dalles 
is inconsistent. 
solution We find that 
A 9'=3 
AI=|2 3 1] =0, JAl= 
Dy Or 2 


Show that the s 
System of e i 
quations 


YN BD 


NNO 


9 

3 ‘ 4 6 3 
1) =0, |AJ=|2 7 1 

~ De] ee. 


Matrizes tihunin 
atrices and Eigenvalue Problems 3.15 


6. 


since |A| = 0 and |A,| # 0, the system of equations is inconsistent 


gxample 3.7 Solve the homogeneous system of equations 


me 23s: 
2 3 —2 y 
A 7 Ale 


Solution We find that |A| = 0. Hence, the given system has infinite number of solutions. Solving 


the first two equations 


bal 


H 


we obtain x = 13z, y = — 82 where z is arbitrary. This solution satisfies the third equation. 
Exercise 3.1 
[ea 2ead! SD 
1. Given the matrices A=|-l 1 1), B=|3 0 -l\, verify that 
Oe Te e8 ot ll 
(ii) |A + Bl #\A\+ [BI 


(i) |AB| =|Al BI. 
2. if AT=[1, —5, 7], B 


, 


3. Show that the matrix A = 


4. Show that the matrix A = 


find A’. 


T 
(3, 1, 21, verify that (AB) = BA- 


joa *2 
2 satisfies 
5: ee 
| 


-3) satisfi 
3 


ip 


2 : = 
the matrix equation A? —4A=S5I= 0. Hence, find A - 


i r : : 4 J] ] = 0. Hence, 
es the matrix equation A> -6A + 5A 
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13. 


14, 


15, 
16. 


17. 


Il 
. For the matrix A=|-2 3 I}, 
aS 


0 1 -l 

4 t 

_ For the matrix A = | 4 3 4/|, verify tha 
3 3 4 


ts 5 =I ; =] 
(i) faqj (A)! = adj (A"), (ii) [adj (A) = adj (A). 


| =2 
verify that 


() (A= (A), (i) (Ay =A. 


For the matrices A =|-l 2 


| Se aa 
1 3 4], verify that 
t= 2 20 9 


(i) adj (AB) =adj(A)adj(B), (ii) (A+B) |#A'+B". 
For any non-singular matrix A = (a,) of order n, show that 
(i) |adj(A)|=|A]"", (ii) adj (adj(A))=|A|"* A. 


For any non-singular matrix A, show that |A =) |Al. 


. For any symmetric matrix A, show that BAB’ is symmetric, where B is any matrix for which the 


product matrix BAB’ is defined. 


. If A is a symmetric matrix, prove that (BA!) (A |B’)! = I where B is any matrix for which the 


product matrices are defined. 


. If A and B are symmetric matrices, then prove that 


(i) A+Bissymmetric,  . (ii) AA” and A7A are both symmetric, 
(ili) AB— BA is skew-symmetric. 
If A and ¥: are non-singular commutative and symmetric matrices, then prove that 
(i) ABT, (ii) AB, <i) Avtp"! 
are symmetric. 
Let A be a non-singular matrix. Show that 
. . 2 
(i) if I+A+A + + A™=0, “then 4= =A" 
i) af 4 2 
MOT Aral. + nearer, Pian ating qin 


Let P,Q and A F a 
be non-singular Square matrices of order n and PAQ = I, then show that 4 


Ifl—A j i 

Aisa non-singular matrix, then show that 

| GA) Sts As 42 
assuming that the series on the right ha ms 
For any three non-singular matrices 


nd side converges. ip ta 


A, B, C, each of order n, show that (ABC) | = € 
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solve the following system of equations: Sand Eigenvalue Problems 3.17 
sc Ss ‘ 


gex-y* z=2, x+3y-zH5 Ee 
i" OF ee a 2x + 4y 2 = Hf eee ' 

20. _xt+y + 2z=2, Bt = eR cae ee 
mo axel, Set peg, ye cee 3y +42 = 6. 


oD: Determine the values of k for wh; 
which the s 
System of equ 


6 fag ations 
cy Zz 0, kx + 3y ee, kz = 0 


has (1) only trivial solution, (ii) non-trivial sol sO a 
: : solution. 
23. Find the value of @ for which the System of equati 
2(sin 0) x +y—-22= SS 
( ) & 2z=0, 3x + 2(cos 26)y +32 =0 
has a non-trivial solution. Bi ee ee ace 


24. If the system of equations x + ay +az= 


0 ES 
non-zero and non-unity, » bx+y+bz=0, cx +cy+z=0, where a, b, c are 


has a non-trivial solution, ‘then show that 


a h c 


25. Find the values of A and g for which the system of equations 
x+2y+z=6, x+4y+3z=10, Kay AP 


has (i) a unique solution, (ii) infinite number of solutions, (iii) no solution 


Find the rank of the matrix A, where A is given by 


a gle PaSDY tg 
34 dauiy om, 
26. ; DT soe) 28 . 
Cee | aihignodigitiie 
> 6 8 
samy 5 
0 c, —-b a 
i ond © q% 5 YS ; is 
29, hip: og var si vo BO.AGa)i| LD Eat Bh nD yess et os ; fal alcibe cot Oy ieeily 2 
F beget eg eel ee ted 2 


‘tian matrix, then iA is a skew-Hermitian matrix and if A is a 
tian matrix. 
> ©, then I + A is invertible. 


31. Prove that if A is an Herm 
skew-Hermitian matrix, then 7A is a Hermi 


32. Prove that if A is a real matrix and A" > 0 as n 


33. Let A, B be » Xn real matrices Then, show that 
(A) + B Trace (B) for any scalars o-and B. 


(iii) AB — BA =1 is never true. 
and ‘C? ='0, then show that 


(i) Trace (@A + BB) = & Trace 
(ii) Trace (AB) = Trace (BA), 
A= B + Cc, BC = CB 


e integer p- 
such that 4); = d,i#j and a; 


34. If B, C are n x n matrices, 
A’*' = B'TB+ (pt 1) C) for any positiv 

35. Let A = (a,) be a square matrix of order ”, 
(Al =e =a)" fo +f - 1)d). 


=¢, i=. Then, show that 
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ermitian, skew-Her mitian or no 5 


' ge . H 
Identity the following matrices as symmetric, skew-symmetric, F 


1 2 3 abe Ore 
36 2 je! 37.16 d e}. 38. | -6 0 e 
aouear 6 a eof e 0 
l 2+41 1-i I 2+4i I1-i 0 2+4i {-; 
BO. 12—47  =5 F=sr/040.)|-24 41 © =5e 325i]! 41./-2+46 0 3-5) 
1+f 345i 6 -l-i 3-Si 6 =[=7 '-34 57-0 
a ee! 0 -i +i | | j 
Aza lds (ONL 43 oe —2i 0 44, |-] 0 l-i 
Lqtdn0 a a a i sho ltt by 2 
1 2 -7 
45. |-2i i 1 
i gage’) 
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14 Solution of General linear System of Equations 


in section 3.2.5, we have discussed the matrix method and the Cramer’s rule for solving a system 
af n equations in n unknowns, Ax = b. We assumed that the coefficient matrix A is non-singular, that 
is |A| # 0, or the rank of the matrix A is m. The matrix method requires evaluation of n> determinants 
each of order (n — 1), to generate the cofactor matrix, and one determinant of order 7, whereas the 
Cramer’s rule requires evaluation of (n + 1) determinants each of order n. Since the evaluation of 
high order determinants is very time consuming, these methods are not used for large values of n, 
say n > 4. In this section, we discuss a method for solving a general system of m equations in n 


unknowns, given by 


Ax=b (3.28) 
a, 2 Uy by, xy 
a ax aca. Gs) b, x 

where Avs sacle 3 27) pel Alexa l 2 
Any) An2 <: inn OF xX), 


aré respectively called the cvefficient matrix, right hand side column vector and the solution vector. 
The order of the matrices A, b, x are respectively m x n, m x | and n x 1. 


The matrix 
Qa, 2 -+ An b, 
Gai, atin ac wear i> 
(A| by=] alte (3.29) 
an G2 °° Gin by 
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is called the augmented matrix and has m rows and (n + 1) columns. The augmented Matrix q 
aoe 


completely the system of equations. The solution vector of the system of equations (3 28 Tibes 
x,,) that satisfies all the equations. There are three possibilities: <8) ig a 


n-tuple (X), X, ---, - 
(i) the system has a unique solution, 
(ii) the system has no solution, 
(iii) the system has infinite number of solutions. 
The system of equations is said to be consistent, if it has atleast one solution and consisteny if: 
has no solution. Using the concepts of ranks and vector spaces, we now obtain the necessary, sf 
sufficient conditions for the existence and uniqueness of the solution of the linear System ‘ 


equations. 


3.4.1 Existence and Uniqueness of the Solution 


Let V,, be a vector space consisting of n-tuples in IR"(or C"). The row vectors R,, ee R,, of the 
m X n matrix A are n-tuples which belong to V,,. Let S be the subspace of V,, generated by the rows 
of A. Then, S is called the row-space of the matrix A and its dimension is called the row-rank of A 


and is denoted by rr(A). Therefore, . 

row-rank of A= rr(A) = dim (S). - (3.30) 
Similarly, we define the column-space of A and the column-rank of A denoted by cr(A). 
Since the row-space of m x n matrix A is generated by m row vectors of A, we have dim (S) <m 
and since S is a subspace of V,, we have dim (S) <n. Therefore, we have 

rr (A) S min (m,n) and similarly cr(A) < min (m, n). (3.31) 

Theorem 3.6 Let A= (a,,) be an m X n matrix. Then the row-rank and column-rank of A are same. 
Now, we state an important result which is known as the fundamental theorem of linear algebra. 
Theorem 3.7 The non-homogeneous system of equations Ax = b, where A is an m Xn matrix, has 
a solution if and only if the matrix A and the augmented matrix (A | b) have the same rank. 


In section 3.2.3, we defined the rank of m x n matrix A in terms of the determinants of the 
submatrices of A. An m X n matrix has rank r if it has at least one square submatrix of order r which 
1S non-singular and all square submatrices of order greater than r are singular. This approach is very 


time consuming when 1 is large. Now, we discuss an alternative procedure to obtain the rank of a 
matrix. 


3.4.2 Elementary Row and Column Operation 
Th ’ ; 
e following three Operations on a matrix A are called the elementary row operations: 


(i) Interchange of any two rows (written as R. ~ R) 
i Tae! 


(ii) Multiplication/division of any row by a non-zero scalar (written as O&R;) 


(ili) Adding/ subtracting + ak; 


ae - a scalar multiple of any row to another row (written as R; — fi ents 
@ multiples of the elements of the jth row are added to the corresponding elem 


of the ith row. The elen ts 0 
; : nents of the jth row : eas, the elemen 
ith row get changed), J remain unchanged, whereas, 
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© not chan 
s said'to Ge is the Tow-rank of A as they do not change 
Ai shite ve equivalent to a matrix A if the matrix B 
uence 
erve that of elementary row Operations. Then, 
(i) every matrix 18 TOW equivalent to itself 


(ii) if A is row equivalent to B, then B is row equival 
alent to A. 


is operations change the form of 4 Be 
qhes the row-space of A. A matrix B ; 
can be obtained from the matrix 

we usually write B = A. We obg 


«iy if A is row equivalent to 
(iii) B and B is row equivalent to C, then A is row equivalent to C 


The above operations performed on columns (that is column i 
ni 


column operat ions. n place of row) are called elementary 


343 Echelon Form of a Matrix 


an m Xn matrix is called a row echelon matrix or in row echelon form if the number of zeros 


preceeding ae non-zero entry of a Tow increases tow by row until a row having all zero entries 
(or no other ¢ DS OM is possible) is obtained. Therefore, a matrix is in row echelon form if the 
following are satisfied. 

(i) If the ith row contains all zeros, it is true for all subsequent rows. 


(ii) Ifa column contains a non-zero entry of any row, then every subsequent entry in this column 
is zero, that is, if the ith and (i + 1)th rows are both non-zero rows, then the initial non-zero 
entry of the (i + 1)th row appears in a later column than that of the ith row. 


(iii) Rows containing all zeros occur only after all non-zero rows. 


For example, the following matrices are in row echelon form. 


b 2aa 4 
1 Sdwi7 biol BS "9 ik 
0 5 4 14,40 0 3 5419 9 0 of 
0 00 9} 10 00 8), 9 9 9 


trix. Assume that a), # 0. If a,;, = 0, we interchange the first row 


i iti -zero. Using elementary row 
t in the (1, 1) position as non 
with some other row to make the elemen aa 
Operati duce the matrix A to its row echelon form (elements of first column be 5 
lons, we reduc ? 


are made zero and so on). 
made zero, then elements in the second column below 4x 


i matrix. 
Similarly, we define the column echelon form of a 


Let A = (a) be a given m Xm ma 


trix A gives the rank of 
: row echelon form of a ma 
Rank f non-zero rows In the — and the set of the non-zero 
the al ne ee es of the row-space Z —— 
: ; row-space- 
tows in the row echelon form gives the basis of the 


a ie 
Similar results hold for column echelon matric 


Remark 13 


() If A isa square matrix, ; ar matrix. 
column echelon form is 2 lowe! triangul 
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(ii) This approach can be used to examine whether a given set of vectors are linearly ind 
or not. We form the matrix with each vector as its row (or column) and reduce it toe 
(column) echelon form. The given vectors are linearly independent, if the TOW eche| the Tow 
has no row with all its elements as zeros. The number of non-zero rows js the ahinene fo 
the given set of vectors and the set of vectors consisting of the non-zero rows i. Tow 
Sis, 


Example 3.25 Reduce the following matrices to row echelon form and find their ranks 


lo 3 al 
to 3 § 
: a Dor ihyt4y 25 
(| eee et (ii) igetSe wha, | 
—2 8 2 
8 1 14 17 


Solution Let the given matrix be denoted by A. We have 


1 5 Poa 
(Gi) Aa al ea ; : sea tok maa . 
. [Posh ew IRS BA el eat a RM a 
Se 0 14 1D 0 0 © 


This is the row echelon form of A. Since the number of non-zero rows in the row echelon 


form is 2, we get rank (A) = 2. 


i: oars se pall ica i 7 feoe2 or 
ean ee SIE (OES MI RSS) BREE E [OTB tas 
Oe i Sa IN) SD NR oe ee 

Seale a4 maliall OQ iF jl) 15 0 ®@ oO @ 


Since the number of non-zero rows in the echelon form of A is 2, we get rank (A) = 2. 


Example 3.26 Reduce the following matrices to column echelon form and find their ranks. 


Battal wii IP eee rey 
oF HEL ag: Sah  I-l =3i=3 
Oy, ae zae (1) ibadatort) 2 

A IS thir =i hel ee esee 


Solution Let the given matrix be denoted by A. We have 


[3 ii 3 0 0 3 0 0 
Pome = (3-0 
(i) A= Bo Cea 5/3 “513 aeeete 5 

er aie Gale3. Sem Aiwa Biss] DrFan peer shar KANTB 0 
a1 is ales Ves (5) pee 
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the column echelon for 6 
since the c0 orm of A has two 
non-zero columns, rank 
s, rank (A) = 2. 

| 1 -l 

eqe MOS CoG | : Deea0 ft 00 
ay A= | 1 0 sel ast eo te SAE) Cet 2G, = 2 0 

[* =1 3 Car “E02 1) Cyt Cy Whit Shed” OT 

i pes oa! 1 20 0 


Since the column echelon form of A has 2 non-zero colum k (A) =2 
: ns, ran = Hs 


Exami I r ; : 
at a whether the following set of vectors is linearly independent. Find the 


example : : 
nd the basis of the given set of vectors 


imension 4 
f 0, 234s 2,0 1, ~ 2), B24, 
Gi. G, {, 0, 1), 4 Ls €1,:1, 1,5 6.0, 1); 
ia 2 3, 6, -3, 4), (4, 2, 12, —3, 6), (4, 10, 12, —9, 10). 
solution Let each given vector represent a row of a matrix A. We reduce A to row echelon form. 
ifall the rows of the echelon form have some non-zero elements, then the given set of vectors are 


linearly independent. 


(23 Apel ji “page A yarn arse: 
; ; > 2k, 
quae |2. 0 452 pare |0 Tae sme10 CRs eee = 
3.9 A he ae. (OS G20. .0 0 


of A are not non-zero, the given set of vectors 
ws is 2, the dimension of the given 
the set of vectors {(1 2 3 4), 


Since all the rows in the row echelon form 
are linearly dependent. Since the number of non-zero TO 
set of vectors ig 2. The basis can be taken as 


(0) 4,55 Be 


fies ie Oa R HrioMe 10 1 pencanon) 
es (plana alge 
ie QrenO wa Oe ae 
(i1) K= a4 1 1 1 R,+R, a 0 2 1 2 R, R,; 0 0 ; 0 A a 
er ee eet a) 0 
od oo) A 4, Og! 
il 2 
alin? Sila jo 2 } 
A Ncdierd: Ra hs lo) Opalae” 
(Om eae 
Ondii2rwA 0 
i t of vectors are 
3 f A are non-zero, the given se 
anes ane ot at given set of vectors 1S 4. The set of vectors 
onan : he given set itself forms the 


0, 0, 0, Dt or t 


> 


linearly independent and 
A(l, 1, 0, 1), ©, 2 
basis, 


” 


2), (0, 0, ’ 
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cme o 3 ROR, 
(iii) A= de 2 12 -3 6 R,-2R 
AeO) 1 Zep oumlO 


Since all the rows in the ec 
linearly dependent. Since the nu 
of vectors is 2 and its basis can 


3.4.4 Gauss Elimination 


Consider a non-homogeneous s 


Ax=b 
ay, 442 
a a 2 
where A= 2! : 
Amy m2 


We assume that at least one element of b is 
m X (n + 1) as 


ay 12 
44, 42 
(A | b) = 
am An2 


=10 


helon form of 


mber of no 
be taken as the set (2, 3, 6, —3,. 4); ©, — 4, 0,3 4) 


not zero. 


Qy, 


Ax 


a 


mn 


2 3 6 =3 4 Be 
=f AOS 1 SOAR Rae qO edie a 
094 s0m 3" S92 On soma | 
A are not non-zero, the given set of Jat 
n-zero rows Is 2, the dimension of the te ate 

Ven Set 


Method for Non-homogeneous Systems 


ystem of m equations in 7” unknowns 


(3.32) 
by x 
b>. 
,p=|2)x=] 
by Xn 


We write the augmented matrix of order 


by 
by 


b 


m 


and reduce it to the row echelon form by using elementary row operations. We need a maximum of 
(m — 1) stages of eliminations to reduce the given augmented matrix to the equivalent row echelon 
form. This process may terminate at an earlier stage. We then have an equivalent system of the form 


a; 42 


0 a) 


(A|b)=] 0 0 


WwW. 
here r < m and a, # 0, a, #0 
2 Sera 


Qy b, 7 

Ap, b, 

am | BF iy 
0] Bei 
0| 5b, 


- > 4, * 0 are called pivots. We have the following cases 
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Then. rank (A) # rank (A | b ‘ 
ermal and r=n ) and the s rity Orso, say by 
; (the number uati "> Om ate Not Zero 
n this case, rank (A) = rank (A | ne Columns j ‘ations has no soluti | 
. a in A * won. 
we solve the nth equatio )=n ) etl : 
n for x and the r+) b. 
is called the back substituti X,, the (n — | § 42> «-+> O,, are all zeros. 
ution method )th equation f 
f Or x,_; and so on. This procedur 
e 


(a) 


ystem of eq 


€ ions 


For example, if we hav 
ave 10 equati 
“der 10 X 6. When r enone 
0 hen rank (A) = tank * In 5 variables. th 
(A |b) = 5 te i tae es augmented matrix is of 
as a unique solution 


* 
Let 1 <n and Dieys Dos b, 
49) +». D. are 
m are all ze 
ros. In this case, + unkn 
} owns, x), x 
lin 229 325 x;. can be 


determined in terms 
of so 
Bete fotcza lee B eee (n =r) unkn 
aaene owns 
(n — r) parameter family of ae for x, and ee i an SMR RC He 
solutions, that is j on. In this 
) 18 infinitel case, We obtain 
y many soluti a 
ons. 


remark 14 
(a) We do not, normally us 
: e column 
equations. iinet We Kee elementary operations in solvi i 

af equations mt er o columns, the order of the oes the linear system of 

‘Keeping ltrack at the sider Serum nowns in the given system 
nknowns is quite diffic 
ult. 


Gauss elimination method may be writte 
n as 


(b) 
(Ae 
Tow operations lc). 
The matrix B is th 
Be, ie Ee echelon form of the matrix A and c is the new right h i 
e solution vector (if it exists) using the back Wena cies 
ion method. 


T 


Gauss elimination method can be used to solve p systems of the f 

Ax=b, which have the same coefficient matrix b i . a a Be ie 

ae ix but different right hand side column v 

| e augment ed matrix as (A |b, Dosee> b,)> which has m rows and (1 + het 

Using the elementary row operations, We obtain the row equivalent s 3 m Apres 

+13) Cp)» where B is the Tow echelon form of A. Now, we solve the ene aie ee 
1? 


B = oo ; 1 1 
X=, ..., BX = p> using the back substitution method. 


a. 
ee 


change this row 


Remark 15 
t. We normally 


0, then we inter 
pivot elemen 
itude. 

ble to compute the 
umber is called 
tion method for 
oximately 


ement becomes zer 
btain a non-zero 
omes largest jn magn 
on. It is possi 
divisions. This n 
he Gauss elimina 
tion count 1s appr 


ation, the pivot el 


a 
(a) If at any stage of elimin 
with any other row below as 


a choose the row such that the piv 
For ann Xx n system, we Ted 
total number of additions, SU 
the operation count of th 
aes an n x n system 18 7 


wire ( s of eliminati 
btractions. ‘plications and 
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Example 3.28 Solve the following systems of equations (if possible) using Gauss elim 


HI method. 
i De 1) [x 4 2 Cau sg 3 
alt 21 2|tyl = |-2) igh al Decl Lil et | dul, 
aie SUE 2 Anoreas See Joel 3 
| easzail gett thatoa 
(cid) | SU dP 
3 SAS 


Solution We write the augmented matrix and reduce it to row echelon form by applying elementary 


row operations, 


| Peis HA lieageas bell Lani 

; (i) (A]b)=| 1-1 2]-2 is af im —3/2  5/2|-4| R, + 5R,/3 
| awn 4\|\ 72 at Qmehs/2.-/22a72) 4 

| 2 foe al 4 


0 -3/2 5/2| —4/}. 
0) 0 8/3)-8/3 


uv 


Using the back substitution method, we obtain the solution as 


z=-—-7, o z=-l, 


2x+y-z=4, or x=1. 


tm 
il 
le 


Therefore, the system of equations has the unique solution x = 1, y = 


2. 0 a3 oa ae pment with’ 
ij R-R/2 = 
(li) (A|b)= 15 8 R - =~ a | 179) 419 R; -2R>= a oi 1/2 — {/2!- 

AuiriacSi ged wiles PO nie Beko bbe ad a ear 

We find that = t 0 
solution. at rank (A) = 2 and rank (A | b) = 3. Therefore, the system of equations has 
ares taal | 
= Pa a 1} J 
(ii) (A|b)=]2 1 ~yJo | -2% Lu Sead 3 
pomp S19 bo 3e-310l B,— Ra lO.. Sorry 
5 —? 2 5 3 SR, 3 2 0 
aha —3/0 0 0 9 
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: i 


Th e 
ve Pe : . numb 
. eatisfied for all values of x, y. > er of soluti 2 
4 ia the first equation, we det From the second Seis We find that the last equation 
Fr X-~ytz uation, we get 3y — 3z = 0, or y = z. 


=|, y= 2 and z is arbitra = ore ak 
: . ry. Therefore, we obtain the solution 


29 Solve the followi 
ample 3 ms System of equations using G 
5 a seet a ae : 
) Ape es 9z+6w =0 (i koe ; uss elimination method. 
ax + 3y + 32 + Ow = 6 Yn 5 


2 syria 
4x —2\y —39z — 6w = —-24, X-3y+z=0 


Sx +y=57 =4 
3x + l4y ~ 122 = 5. 


golution We have 


Sota ways Boe af 6 
eee 0 


4 
(i) (A | b) 2 3 33 6, aeRO 
a) og ade otaeie SAS 9/2 15/2 3) 6|R,+4R, 
0 -18 -30 -12|=24 ; 
A =3) 99) 961/10 
=110) OY 2 aIS 2S nO) 
0 0 0 0) 0 


The system of equations is consistent and has infinite number of solutions. Choose w as 
arbitrary. From the second equation, we obtain 


IS Ba: 15 1 
me 5 y Al Ww = 2| ae 5z —2w). 


From the first equation, we obtain 

4x = 3y + 9z = 6w=4 — Sz 2w + 92 ow = 4 47 — Bw 
or X% al tezae Ze 
Thus, we obtain a two parameter family of solutions 


x= l+2-2w and y= (4-52 —2w)/3 


where z and w are arbitrary. 
Vea paar ne 2 


je ee al 
| Ree i eee | Ue ae 
(i) (A |b) = A esere (1 R,-5Ri ~\9 -9 5|-4| Rat 8Ro!7 
y R, -3R ealieig 
gen jap eee page 
2 | ae : 
i= 5 —A lp, =R/5-\9 o 10/7 —10/7 
- _ 10/7 
0 10/7 | —10/ ast 0 0 
0 ai7| -2/7 


£ ——s——| 
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ed for all values of x, y, Zz. From the third equation ay. 
1. Hence, the system of equations has a Mei 
(24R, — 7R2 + R3)/5, the last equation jg aa Solut 
Unden 
t 
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The last equation is satisfi 


Back substitution gives y = 1, ¥ = 


x=1,y =I andz=1. Since Ry = < 


3.4.5 Gauss-Jordan Method 


In this method, we perform elementary row transformations on the augmented matrix [A} 5 
| and 


reduce it to the form 


Elementary 


[A | b) ————— [I | ¢] 


row operations 


where I is the identity matrix and ¢ is the solution vector. This reduction is equivalent to findj 

solution as x = A’'b. The first step is same as in the Gauss elimination method. From sea in 
onwards, we make elements below and above the pivot as zeros, using elementa ue | 
transformations. Finally, we divide each row by its pivot to obtain the form [I | c]. alteat 4 
every step, the pivot can be made as | before elimination. Then, ¢ is the solution vector, ee 


This method is more expensive (larger operation count) than the Gauss elimination. Hence We do 
not normally use the Gauss-Jordan method for finding the solution of a system. However. chien 
is very useful for finding the inverse (A!) of a matrix A. We consider the augmented matrix [A | I] 


and reduce it to the form 


Elementary 
A. ———— ee sa ned 


row operations 


using elemeniany, row transformations. If we are solving the system of equations (3.28), then we 
have x = Ab, and the matrix multiplication on the right hand side gives the solution vector. 


Remark 16 
If any pivot element at any stage of elimination becomes zero, then we interchange rows as in the 


_ Gauss elimination method. 


xample 3.30 Using the Gauss-Jordan method, solve the system of equations A x = D where 


0 
A=|2 1 -3}, and b=1|4}. 
] l | | 


Solution We perform elementa 


ro : . the 
form [I | c]. We get ty row transformations on the augmented matrix and reduce I! th 


] —| ] 0 
we GAT 3 oul 
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| 


10 <A SS7aianelie Rae le eon 
2 TA ES ORS 0 1 ~5/3| 4/3) R,/(10/3) 
0 0 10/3|-5/3 


S 
r=) 
ee] 
Nm 
+ 
Nn 
Re) 
a) 


Hence» the solution vector is 
x=[l 1/2 1/2)" 


F =I i 2 
pxample 3.31 Using Gauss-Jordan method, find the inverse of the matix A=| 3 —1 1 
==) ees) 4 


Solution We have 


=] .1.2,) 1.0.0 
(A|T)=| 3 -1 11/0 1 0}. 
—| 3 4)0,0 1 


The pivot element a,, is —!. We make it 1 by multiplying the first row by —1. Therefore, 


reser alls Oa | - =2\-1-0 0 
Ry -3R 
ase 3 Stee 0 Olek chew dO7e? 7\ 3 1 0\R/2 
oy ail 0) Oleg ie eel 
3/2| 1/2: 1/2 0 
ay 2 2)| lee i 
Li bn bi 6 2 Oe 1/2 0\ CR)/5 
~ 10. Ip Wale 4 RIP ire ee ae 
oil: Siena 0 0 =5 


0 fo of -7/10 | 2/10 3/10 
0 3/2|1/2 1/2 R,-3R3/2 2} 9 1 0)-13/10 —2/10 7/10). 
00 | 


— 


wo des T/2SIe. ee 0} a, —7Rs/2 his ee ills 
5 -l 
0 One alliaie 1/ 
ag. Digee 
1 Eh did 
t= — —13 rd 
Hence, A | 3 2 4 
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3.4.6 Homogeneous System of Linear Equations 


Consider the homogeneous system of equations 
Ax=0 
where A is an m X n matrix. The homogenous system is always consistent since x 
solution) is always a solution. In this case, rank (A) = rank (A | 0). Therefore, for the eve (trivig 
92 io 9 . Ee Se. g 

system to have a non-trivial solution, we require that rank (A) < n. If rank (A) =; < i He 
. . . . . 4 > Q . 
an (7 — r) parameter family of solutions which form a vector space of dimension (n — ras ( Dain 
n~ 


parameters can be chosen arbitrarily. 


The solution space of the homogeneous system is called the null space and its dimension is cal 
the nullity of A. Therefore, we obtain the result a 
rank (A) + nullity (A) = ” (see Theorem 3.5). 


Remark 17 


(a) If x, and x, are two solutions of a linear homogeneous system, then ax, + Bx, is also a 
solution of the homogenous system for any scalars a, B. This result does not hold for 
non-homogenous systems. 

(b) A homogeneous system of m equations in n unknowns and m < n, always Possesses 
non-trivial solution. 


Theorem 3.8 If a non-homogeneous system of linear equations Ax = b has solutions, then all these 
solutions are of the form x = x) + x, where Xo is any fixed solution of Ax = b and X;, 1S any solution 


of the corresponding homogeneous system. 
Proof Let x be any solution and x, be any fixed solution of Ax = b. Therefore, we have 
3 Ax=b and Ax) =b. 
Subtracting, we get 
Ax — Axy=0, or A(x — x9) = 0. 


Thus, the difference x — x9 between any solution x of Ax = b and any fixed solution %) of 
Ax = b is a solution of the homogeneous system Ax = 0, say x,. Hence, the result. 


Remark 18 


If the non-homogeneous system Ax = b where A is an m x n matrix (m > n) has a unique solution 
a is rank (A) =n, then the corresponding homogeneous system Ax = 0 has only the trivial solution, 
that is x, = 0, 


Example 3.32 Solve the following homogeneous system of equation Ax = 0, where A 1s given b) 


(2 | limes ae) el 
(i) —1], GO =e Gite eee te 4, 
a el 2 Bee 6: I 


Find the rank (A) and nullity (A). 
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ite | Cc ente oy es 
tio” We write the augmented matrix (A | Q) and reduce j 
$0 €duce it to row echelon form. 
2 1| 0 > 
() Al Ry —3R /2 ~|9 -3/2|0| R 
g (3 oN /2 ; 3+R)/3 =|0 —3/2/0]. 
ie 0 lo 
since, rank (A) = 2 = number of 
Hence, nullity (A) = 0. unknowns, the system has only a trivial solution. 
COs |e 2 SMO Salo 
(1 << R,-R, =1 5 2/0)R,-3R, =|0 —1 © 2) 0). 
=lo Ae 0 =3 1 4io 0 0 -2/0 
since rank (A) = 3 = number of unknowns, the homogeneous system has only a trivial 
solution. Therefore, nullity (A) = 0. 
My ts NS eee i WeSiiy MO Pe Sih HO 
Gi) (AL 0)=|2 3 Ol, a, = (C 1 3 DOR | 1s 210 
eo) 2 9 (0 =i =3 <2) 0 0 0 0 ojo 


Therefore, rank (A) = 2 and the number of unknowns is 4. Hence, we obtain a two parameter 


family of solutions as x2 = —3%3 2x4, Xy = — Xp + Xy — Ny = 44s + x4, Where x5 and x, are 


arbitrary. Therefore, nullity (A) = 2. 


Exercise 3.3 


Using the elementary row operations, determine the ranks of the following matrices. 


De Ouel ee Ye? 
has ie Ol yo ot 2 5 (eles eal. 
M? t 3 6 = MU 3 Fira 2 
eee 
i dep irae Patel 5 ath 
Aiomaet 4% 5 2? 2 Clie, nee a 
Ae e1g 35 A Say gr g 13 14 
yO IO - 
ee ee Te alow Ao 5 
Te foe uaa mbbasiee s Pally see Bee 
7 11 15 19 ogee Pesee 6 
9 15 21 27 eli 
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ees a 
Sue el | 
1. tg A 2 se 
Seuss 5 
Using the elementary column operations, determine the rank of the following matrices. 
l -] 23), gael E2 1 Ds ©) 
AileellnleesO) Sl 12. 3 call Ile 13. |-1 1! 5 
fe Si 3 Sed 4 3 A 8 
(Peo? =3 De 30 0 4 
14 2 ] —-l 15 3 1 2 -l 1 
ek al eet ga 2-221 
5 45> Kea. isl 5 
Determine whether the following set of vectors is linearly independent. Find also its dimension. 
5 KOx 2 es She Dh CAD 


16. {(3, 2, 4), (1, 0, 2), . -1, -D)- 
18. (2, 1, 0), (,-1, 1), 4 1, 2) 2 3, Dt 194 (202) 1), Ont P(e er 
PUG laa, Dau, I —1 2) 

(1, 1, 1, Ys Cl. 1, 1, Ds C1 1, Ds Cy 1, OD): 

(1, 2, 3 1), 2 1, WD), 4, 5, 5, 3), G, 4 1 3) 

(1, 2, 3, 4), (0, 1, -1, 2), (1, 4, 1, 8), (3, 7, 8 14)}- 

(1, 1,0, 0, (1, 1, 1, D, 4 4, 1, D.C, 0 ® DI. 


21. 
22. 
23. 


24. 
HS, 12, 2.0, Dy (Ch ta 25 Wh (3, 0, 4, 0)}. 


Determine which of the following systems are consistent and find all the solutions for the consistent system 


Fee i p) t 4 Flexi ii A Fle: 

Pieiae col. Zula ile if, 2 S89 1S) 2 28) 130 8 aeiiimiaie 
De Sz) D2 BIZ = ee Gililezalate 

cof kl ie Le ed eae die ehOeeaille 
yl=|-4|. 30. {1 2 3/|y/=] 16]. 210 12 allo 
-3 4][z] | 6 1 3 4]|z| {22 1. Hah 
5 314 4 

ie relent ee Se eee ae test NE 
fal 21/71" 10 BS est, Llp ON. BA fie oh 2 gel i 
2. 6a a2 foes iiiet | Bla Teeee elie 

w Ww 
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y|_|0 
Beye Laat Lz tol 
0 


rind all the solutions of the following homogeneous syste 
ms Ax = 


Zee le tne ee = 0, where A is given as the following. 
2 
45 3 Siteu\esl der 3.-11. 5 
1 5 -4 25 ee) ' 38. |4 1 —-10 
yo. 6 
ayy ee Age (alia, Asien rT aS 
ome faire) Ths 41. 
ieeeon 2 1 ape ard aie 
mine ae i of 1 1 
ape e4 pe 3 5 spl nth 
4 0 4 10) 1 43 2 =\| =2 =3 y 1 Tao 8 
eee sy ssy- <3 Tee iOse els cL ates 
mn 2 4 3 5 —| =, =A 5 7 > = 
i ee 
Dail 1 =2 
45. ee Ar 
ee) ee) 


Using the Gauss-Jordan method, find the inverses of the following matrices. 


=e dete 
BBS ih et 
rae 1 3 3 48. a 
AG. lime 3 47 “ate 
Onl 2 i. eat 
| Be 4 
Peele 1 1 
con false 
49, |~ x S01. aaa 
== pe . amon 
fe a eh 
ular. Consider 
- pawn Pe atrix A may be singular oF non-sing 
n. The m 
Let A= (a,) b matrix of order (3.35) 
e a square 
ons 
the homogeneous system of equatl Boat ea AN) x=9 
sig 30 
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where A is a scalar and I is an identity matrix of order n. The homogeneous System of 
(3.35) always has a trivial solution. We need to find values of A for which the homogeneo, tion, 
(3.35) has non-trivial solutions. The values of A for which non-trivial solutions of the ho s 
system (3.35) exist, are called the eigenvalues or the characteristic values of A and the coreg 

non-trivial solution vectors x are called the eigenvectors or the characteristic vectors of A : 8 
a non-trivial solution of the homogeneous system (3.35), then @ x, where @ is any constant oa 
a solution of the homogeneous system. Hence, an eigenvector is unique only upto a constant rai 
The problem of determining the eigenvalues and the corresponding eigenvectors of a square tee 


A is called an eigenvalues problem. 


3.5.1 Eigenvalues and Eigenvectors 


If the homogeneous system (3.35) has a non-trivial solution, then the rank of the coefficient matriy 
(A — AD) is less than n, that is the coefficient matrix must be singular. Therefore, 


aya ap) st Ay 
a. =f} eae 5p 
sean | ri : =(}. (3.36) 
ay| an? a3 One A 


Expanding the determinant given in Eq. (3.36), we obtain a polynomial of degree n in A, which is 
of the form 
P(A) =|A = AL) = C1)" [A" = 0, AT! + cy A? - 2. +:Cl)" c,] = 0 
or M = eA + cn W*— ... + Cl)"e, = 0. G5) 
where C;, C>, ..., ©, can be expressed in terms of the elements a, of the matrix A. This equation is 
called the characteristic equation of the matrix A. The polynomial equation P(A) = 0 has n roots 
which can be real or complex, simple or repeated. The roots A,, A>, ..., A,, of the polynomial equation 
P(A) = 0 are called the eigenvalues. By using the relation between the roots and the coefficients, 
we can write 
A, te A, te. 0 A, = fy a a) + a9 i Any 


AjAn AAs +. +A, Ay = 


AA; ese A, ="C (3.38) 


If we set A = 0 in Eq. (3.36), then we get 
[Matql)te =e = 2a. ei a 


” 


n 


Therefore, we get 


sum i a 
m of eigenvalues = trace (A), and product of eigenvalues = Al. 
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he » see rd . Dect: 
fhe ne spectral radius of A and is agithenk “um Of A and the | 
e 1e 


g), we f pete 
(3.38 ), we find that atleast one of tl by p (A). If |A| =0 Sa eigenvalue in magnitude is 
Ne elgenva » that is the martix is singular then 


es iS ZrO, then |A| = 0, N | 
cig ses VS Note: that 4 ues must be zer 
matrix, then the at if A j PEIORS i 
ular I diagonal elements a IS a diagonal or an u ene ed ifvone:oF the 
s ne matrix A pper triangular or a lowe 
are the eigenv F 
genvalues of A. 


the h 
Omogeneous system (A — A; I)x = 0 for each 
1 


tors. 


4Q 


etermining the eigenvalues A,’s, we so] 
i% SOlve 


_.., n to obtai 
oe btain the Cotresponding ej 
elgenvec 


erties of eigenvalues : 
prop and eigenvectors 
pet A be an eigenvalue of A and x be i : 
results. 'S corresponding eigenvector. Then, we have the follow! 
: , e the following 


_ qa A has eigenvalu 

I g © aA wee eigenvector is x. 

2. A” has eigenvalue A” and the pet 5 (QA)x. | 

Pre-multiplying both sides of Ax = Ax by N wee soa resieh aly HN ae 
Therefore A’ has =e a 7 . HOD 

Pe evel ee aa = pa i ae corresponding eigenvector is X. Pre-multiplying 
A — ki has the eigenvalue A — k, for any scalar k and the corresponding eigenvector is X. 


Ax =ax > Ax-kIx=Ax—kx 


(A —kD)x= A-1x. 


or 
alue 1/A and the corresponding eigenvector is X. 


4, Alf it exists) has the eigenv 
Pre-multiplying both sides of AX = Ax by A ', we get 
aay or A x= (I/A)x. 


A 'AX = 
4 — k) and the correspondin 


Sra = kt)! has the eigenvalue 1/( g eigenvector is x for any 
scalar k. 

6. A and A’ have the same ei 
columns) but different eigenvectors, (s 

7, Fora real matrix A, if at ip is an eigenva 
(since the characteristic equation has real coe 


property does not hold. 
ch gives the relationship of a ma 


genvalues (since a determinant can be expanded by rows of by 
ee Example 3,41). 
Jue, then its conjugate O— 


fficients). When the ma 


jf is also an eigenvalue 
trix A is complex, this 


trix A and its characteristic 


We now present an important result whi 
n characteristic 


“A satisfies its OW 


equation. 
matrix 
Theorem 3.9 (Cayley-Hamilton theorem) Every square 
equation vee pis? ee 
A! At! oe se eee Gal At (-l) Cn Z 
— C 
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Proof The cofactors of the elements of the determinant |A — AT| are Polynomials : 
(n — 1) or less. Therefore, the elements of the adjoint matrix (transpose of the cof wa 


0 

actor m Wepre 

also polynomials in A of degree (m — 1) or less. Hence, we can express the adjoint Atrix) s 
Mat, 


polynomial in 2 whose coefficients B,, B, ..., B, are square matrices of order » having 84 
as functions of the elements of the matrix A. Thus; we can write ng “lemon 
adj(A— Al)=ByAtr + By? +... +B, AB, 

We also have 
(A - AI) adj (A-AD=|A-AIL 
Therefore, we can write for any A 
(A — Al) (B\A"! + BLA"? +... + B, A+ B,) 
=MI-c Art... HEI e AT + Cl)" 1 


Comparing the coefficients of various powers of A, we obtain 


Se 
AB, - B= ¢,1 
AB, — B, = ope | 


ABy | = By= 1)" e,41 
AB,, = (-1)" c, I 
Pre-multiplying these equations by A”, A”~', ..., A, I respectively and adding, we get 
. A”—cA"" +... + Cl)", A + (-1)'c, I = 0 
which proves the theorem. 


Remark 19 
(a) We can use Eq. (3.39) to find A” (if it exists) in terms of the powers of the matrix A. 
Pre-multiplying both sides in Eq. (3.39) by A’, we get 
Av!A? ar Cy Anan! de aa + (1)"¢ AA + (-1)"c,, AUT = AO Ae, 0 
ae Ge! ial aa 
or Ng We [A"> 1 —¢ A"? +... + Gyr! ¢, 21] (49) 


co, 


(b) We can use Eq.(3.39) to obtain A” in terms of lower powers of A as 
ANS GAN! oA"? +. + C1 e I. 
Example 3.33 Verify Cayley-Hamilton theorem for the martrix 
P28 
A=j-1 1 2]. 
Li 4 


| 
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‘) obtain A and A’, (ti) find eigenvalues of A, A? 
; es of those of A, (iit) find the spectral radius of A and verify that eigenvalues of A? are 
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golution The characteristic equation of A is given by 
1A 2 0 
|A—AR=) (2-1 42 s | 
A NSS) = aay a2 Ay} 
oa 
= (1 =A) (A? = 22-3) 2 3) = = + 302 - A+ 3-=0. 
Ihe lo 2 Ol) [ai dea 
Now, SSS 1 DW nt l= 3 4 
TU a) I 6 5 
a 4 ANT i ® | 1 1 
AS = AZAy= || 80) 93 lt) 2) el nO) 
0 Si 1 2 a1 1G Wy 
We have 
= 10 12 1 4 4 12 O 100 
8 3a? — A +3 = a ae 2) islet at A Eeei|G) 
Hh Salli 06 S eS Mie Oo 
OO (3.42) 
0 0 0 


isti nee BoA 3 = 9: 
Hence, A satisfies the characteristic equation — 4 


(i) From Eg. (3.42), we get 


meaty 

lilies ae | 0-0 pe inet 
| 1 2 le = il ae | x3 0 3 
rasan We te 3 PG) \Ome 
From Eq. (3.42), we get 12 

| pol (2 he le 
a&12 12 || Ona pale 8 
sale af ite 1 

A=3a2-A+31=| 0 ? at) oe ; 

0 18 } 
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oots of 


(ii) Eigenvalues of A are the r 
Pore tb Oot Aas bi, 


(Pezey ganas 


The characteristic equation of A’ is given by 


ge 4 
ee RO IG SAT 2A=0 
0. @ dak 


A+ G2= 8-9 =00r + YA=9NAHY=0 


or 
S11, 41 which are the squares of the eigenvalues of A. 


The eigenvalues of A’ are 


(iii) The spectral radius of A is given by 


p (A) = largest eigenvalue in magnitude = max |A,| = 3. 


ile WO 
Example 3.34 IfA=|1 0 1}, then show that A” = A"-2 + A* —I for n = 3. Hence, find 4” 
0 a0 


Solution The characteristic equation of A is given by 


1A.) 0% © 
|A-All= foe l= @ DG adie 0, or AA AFI —0 
Op ALA 


Using Cayley-Hamilton theorem, we get 
Ao-A2-A+1=0, or A’-A’=A-TL. 
Pre-multiplying both sides successively by A, we obtain 
A?-A°=A-I 
A‘ A?=A’-A 


Atl A’? = A™3 ty A+ 
A” At! = A”? ~ A”™3 
Adding these equations, we get 


n ran ee 
A’ — A-= 4! a aa AM Al 2 Ke ee 


° 


Using this equation recursively, we ‘get 
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A" = (An4 we ‘and Eigenvalye Broo aei 
AD) ct 42 
= (A"-6 ee ) Aes I = Ar-4 4 2(A2 
pba? veges eee 
| eA aes (ACen) 
= AC a a 
a n—~ 2) (a2 _ heh 
2 [ieneunar hs 
‘ = 50, we get / > A 5 (n 21. 
gbstitutine n ? g 
U 
-O @ 
A” = 2542 9 10 0] F100 
=r = 
P= 25/1 1 ol-2alo a ol=l2s-4 0 
SOUS NO Oa) as Oa 


gxample 3.35 Find the eigenvalues and the corresponding eigenvectors of the following matrices. 


ans ja 10 0 
OSS Wee ale wa-| (iii) A=|0 2 1). 
2 0 3 
Solution 
(i) The characteristic equation of A is given by 
-|1-A 4 ; 
|A —AI|= 3 22) =0 or A-3A-10 =0, or A=—2, 5. 


Corresponding to the eigenvalue A = —2, we have 


3 4\( x, 0 £ ey 
aernx= [5 eae or 2 eines 0 or x, ae 


Hence, the eigenvector X is given by 


ule be a [ = 
aS xX X i 2 
nt multiple, we can take the eigenvector as [- 4, aikic 


i is uni a consta 
Since an eigenvector is unique upto 


; = 5, we have 
Corresponding to the eigenvalue A=: 


(A-SD*X-| 3 23)\% 


pread yond 
= (X), “2 
Therefore, the eigenvalue is given by x (1 


a of A is given by 


2 = 0, 
1A }-° or Ne, 
Aca _1 1-4 


(i) The characteristic equatio Bets 
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or 


Both the equations reduce to — x) ~ 1X5 


Corresponding to the eigenvalue ,Z=1+ i, we have 


S| 0 
[encler’) I] x= B EAL 


+7 1X, + Xo = 0 and i xy are 1X, = 
= (0, Choosing x2 = 1, we get x) = -i, Thee 
g, the 


0. 


eigenvector is X = [-i, ike 


Corresponding to the eigenvalue A 


or 


Both the equations reduce to —x, + 1X5 


is x = [i, 1] 


= | — 7, we have 


ae ppt 0 
aa -anx=| | NEL 


ix, + x= 0 and— x, + iX, = 0. 
= (), Choosing x, = 1, we get x, = i. Therefore, the eigenvector 


if 


Remark 20 


For a real matrix A, the eigen 


values and the corresponding eigenvectors can be complex. 


(iii) The characteristic equation of A is given by 


(Aeon 2A A OR COS Ge Ae Uns Aad 


2 0 3-A 


Corresponding to the eigenvalue A = 1, we have 


0 0 OF}; x, |} 0 
(A-Dx={]0 1 If}x|]/0] or ie, 
2 0 2i/x,{I0 maa © 


We obtain two equations in three unknowns. One of the variables x,, x5, x3 can be chosen 
arbitrarily. Taking x; = 1, we obtain the eigenvector as [—1, —1 Ie 2 


Corresponding to the eigenvalue 2 = 2, we have 


=1,.0 
(A-21)x=| 0 0 
2 0 


O}] x 0 
] XY) = 0 
Tie 0 


ie 


or x; = 0, x; = 0 and x, arbi 
2 arbitrary. Taking x, = 1, we obtain the eigenvector as [0, 1,0 
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orres = 
C 3, We have 


f2om4 6 
28 we x1 1/0 
(A - 31I)x= ik 
X5 0 or x, =0 
A 1) 1 x, |[0 


Choosing x3; = !, we obtain the eigenvector a 
S 


'e rs [ 0, =|, 1 Pr 
gyample 3.36 Find the eigenvalues and the correspondi 
Fou, Ponding eigenvectors of the following matrices 
; 4 lik .@ 
(GVA = 70P oleae hs OO 
Apiyet (ii) A={0 1 of, Gia = lo Tum0 
0 Out 001 


solution In each of the above problems, we obtai 
f , tain the characteristi ; Bre 
the eigenval = cteristic equation as (1 — A) = 0. 
ee EGR aul aa “.t I, 1, a repeated value. Since a 3 x 3 matrix has : sabi 
ipis"unp » Whether the given matrix has 3 linearly independent eigenvectors, or it has 


lesser number of linearly independent eigenvectors 


Corresponding to the eigenvalue A = 1, we obtain the following eigenvectors 


' Ot OF /@ X, =0 
(i) (A-I)x=|]0 0 1|x=/0) or i) 
0 0 0 0 x, arbitrary. 


Choosing x, = |, we obtain the solution as [1, 0, oy! 


Hence, A has only one independent eigenvector. 


One Oliies 
(i) (A-Dx=|0 0 O|]x,|=|0), o a 
og xs. “2a x , x3 arbitrary. 
0 0 0 


x3 0 
Taking x, = 0, x3 = | ENGL Te ae Wey = 0, we obtain two linearly independent solutions 
T 
x, = [0, 9, i, = 01.0, OF. 


In this case A has two linearly independent eigenvectors. 


0 0 O01" 0 
(iii) (A - I)x= 0 0 0 xX =| 0), 
‘to 0 Of[%s 0 


ary values of all the three variables. Hence, we obtain three 


Thi is sati or arbitr 
s system is satisfied f 5. which can be taken as 


linearly independent eigenvector 
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ies | ‘ 


x [1, 0, 0], 2 [0, > Vs 3 [O, O, 1] 


ults regarding the relationship between the eigenvalues of am 
Atrix 


dependent eigenvectors. 
net eigenvalues are linearly independent. 


2. If A is an eigenvalue of multiplicity m of a square matix A of order n, then the ani 


linearly independent eigenvectors associated with A is given by 
where r=rank (A-AD,1spsm. 


We now state some important res 
and the corresponding linearly in 


1. Eigenvectors corresponding to disti 
er of 


p=n—, 


Remark 21 


In Example 3.35, all the eigenvalues are distinc 
linearly indenpendent. In Example 3.36, the eigenvalue A 


t and therefore, the corresponding eigenvectors fe 
= | is of multiplicity 3. We find tha - 


(i) Example 3.36(i), the rank of the matrix A — I is 2 and we obtain one linearly indepen ia 


eigenvector. 

(ii) Example 3.36(ii), the rank of the matrix A — I is 1 and we obtain two linearly independen; 
eigenvectors. 

(iii) Example 3.36(iii), the rank of the matrix A —I is 0 and we obtain three linearly independen, 
eigenvectors. 
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Matrices af 
Ices and Eigenvalue Problems 3.91 


pxercis® Si 
rity the Cayley-Hamilton theorem for the matrix A. Find Al ip cos 
6 , If it exists, where A is as given in Problems 
j to" 
1 0 -4 yoy 
5 ae) a 
ae Cpe = Peal (ee 7) 
=f 4-3 124 ; 2 
‘au Be hell ool Gea 
Ae alls x veer eal eae ne 3 
it Soe ap 1 tj 


pind all the eigenvalues and the corresponding eigenvectors of the matrices given 
Which of the matrices are diagonalizable? 


Ste West a fered 
TNO: 2a 8. |1 —-1 0). Fy 
x le Poe 0S 2 OI 
im 3 eB he Wg Oita 
1, | as) ii, |} b @ (We OF IIe 
-1 3 =i —i | bE ue 
igal a Ow =2 Gy cil gt 
ito @ = fyeCiaal al 

eae ug 2 se oh tos Be4ORe 

"10 0 0 1 Hl oliver Qt 

000 0 = lea ed J eg) a, 

) O OO G1 0 @ 0 

fi: Dm Brat iy Or O 0 Get oO 0 

fe Mia BY 17.0 0 1 0 0 18.0) OOF Ok 

ae ZS OF yo WY OO 00001 

Beta eal? 10000 Vr DhoO O40 


sal 

i i i he matrix P such that P AP 
5 : 24 are diagonalizable. Find t 
Show that the matrices given 1 Problems 19 to 


‘Sa diagonal matrix. 


pene (ee 
4.927 5 he oO 3h 
19 fina 0 4 ecw 
0 2 0). Pees 
ina 
| eee 
5 -06 —9 1 
ea SE, ae 
ees 
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3.92 Engineering Mathematics 


d the corresponding eigenvectors are as given in Problem 
s 


r 2 
0)”, (-1, 0, 1)’, (A, 0, 1)”. Sto) 


Find the matrix A whose eigenvalues an 


25. Eigenvalues: 2, 2, 4: Eigenvectors: (-2, 1, 
_|, 2; Eigenvectors: (1, 1, 0)/, (1, 0, Lye (,. 1; 1)°. 


26. Eigenvalues: 1, ) 
2, 3; Eigenvectors: (1, 2, 1) ee 4) 5 lb, 45 9) 2 


27. Eigenvalues: I, | 
1, 1; Eigenvectors: (-l, 1, 1 es (lb = bys (GBs aby. 


28. Eigenvalues: |, 
-1, 1; Eigenvectors: (-I, 1, 0) eu; =) (ile ibs by, 


29. Eigenvalues: 0, - 


30. Eigenvalues: 0, 0, 3; Eigenvectors: (il; 4 = es (-2, |, 0)LNG20) 1) 


31. Let a 4 x 4 matrix A have eigenvalues 1, —1, 2, —2, Find the value of the determinant of th 
e 


matrix B= 2A+A!-I. 


32. Let a 3 x 3 matrix A have eigenvalues 1, 2, —1. Find the trace of the matrix B = A ~ ,-! ay 


33. Show that the matrices A and P-'AP have the same eigenvalues. 


34. Let A and B be square matrices of the same order. Then, show that AB and BA have the same 


eigenvalues but different eigenvectors. 
Show that the matrices A-'B and BA™ have the same eigenvlaues but different eigenvectors. 


35. 
Ann Xn matrix A is nilpotent if for some positive integer k, A‘ = 0. Show that all the eigenvalues 


36. 
of a nilpotent matrix are zero. 
If A is an 7 Xn diagonalizable matrix and A’ =A, then show that each eigenvalue of A is 0 or]. 


37. 
: ah f : ; 
38. Show that the matrix A = ie a # b, is transformed to a diagonal matrix D = P ‘AP, where 
h 
cos@ —sin@ 2h 
P is of the form P=] . and tan 20 = , 
sin@O  cos@ a-—b 


Let A be similar to B. Then show that (i) A! is similar to B |, (ii) A” is similar to B” for any 
positive integer m, (iii) |A| = |B]. 

Let A and B be symmetric matrices of the same order. Then, show that AB is symmetric if and 
only if AB = BA. : 

41. For any square matrix A, show that A’A is symmetric. 
42. Let A be a non-singular matrix. show that A’‘A | is symmetric if and only if [= (A’). 

43. If A is a symmetric matrix and P-'AP = D, then show that P is an orthogonal matrix. 

44. Show that the product of two orthogonal matrices of the same order is also an orthogonal mattix. 


39. 


40. 


hom ny 


45. Find the conditions that a matrix A = lb my ny| is orthogonal, 


hk m3 nz 


46. If A is an orthogonal matrix, show that |A| ah IK. 


47. Pr ; l 
ove that the eigenvectors of a symmetric matrix corresponding to distinct eigenvalues are orthogom 


48. A matrix A is call fetes =t an and 
called a normal matrix ifAA! = RULAG Ghowiihatthouclemnidant skew-Hermitian 


unitary matrices are normal. 
49, i i 
2 2 a sere. can he diagonalized using an orthogonal matrix, then show that A is symmeltie- 
- Suppose that a matrix A is both unitary and Hermitian. Then. show that A = I 
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